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ABSTRACT 

We report the Chandra discovery of an X-ray cluster at redshift z ~ 1.063 associated with the Compact 
Steep Spectrum radio loud quasar 3C 186 (Q0740-I-380). Diffuse X-ray emission is detected out to 
~ 120 kpc from the quasar and contains 741±40 net counts. The X-ray spectrum of the extended 
emission shows strong Fe-hne emission (EW=412eV) at the quasar redshift and confirms the thermal 
nature of this diffuse component. We measure a cluster temperature of 5.2^|!) g keV and an X-ray 
luminosity L(q 5_2keV) ~ 6 x 10*^ erg sec~^, which are in agreement with the luminosity-temperature 
relation for high-redshift clusters. This is the first detection of a bright X-ray cluster around a luminous 
(Lboi ~10*'' erg sec~^) CSS quasar at high redshift and only the fifth z > 1 X-ray cluster detected. We 
find that the CSS radio source is highly overpressured with respect to the thermal cluster medium by 
about 3 orders of magnitude. This provides direct observational evidence that the radio source is not 
thermally confined as posited in the "frustrated" scenario for CSS sources. Instead, it appears that the 
radio source may be young and we are observing it at an early stage of its evolution. In that case the 
radio source could supply the energy into the cluster and potentially prevent its cooling. 

Subject headings: quasars: individual (3C 186) - X-rays: galaxies: clusters 



1. INTRODUCTION 

Powerful radio sources, that are compact on galaxy 
scales, the Giga-Hertz Peaked Spectrum (GPS) and Com- 
pact Steep Spectrum (CSS) sources, comprise a signifi- 
cant fraction of the bright radio source population (10- 
20%, O'Dea 1998). They are strong candidates for being 
the progenitors of large-scale radio sources (e.g. Fanti et 
al. 1995, O'Dea & Baum 1997, O'Dea 1998), but this 
connection has not been firmly established. Their ra- 
dio morphologies show compact emission on arcsec (VLA 
resolution) scales while on milliarcsec scales (VLBI) the 
sources look remarkably like scaled down large radio galax- 
ies, where the entire radio structure (1-10 kpc) is en- 
closed within the host galaxy. Since the first GPS sam- 
ples have been constructed there has been a clear contro- 
versy regarding their nature (see O'Dea 1998 and refer- 
ences therein). In the evolution model the source size and 
the characteristic spectral break at GHz radio frequencies 
could be an indication of young age, while in the other 
model the radio jet could be frustrated by a dense con- 
fining medium. Recent observations (e.g. measured ex- 
pansion timescales of < 1000 years, Owsianik et al. 1998, 
Polatidis & Conway 2003) give more weight to the evo- 
lution model (Readhead et al. 1996, Snellen et al. 2000, 
Alexander 2000), although there has been no definite ob- 
servational evidence to rule out either of the models and 
both interpretations are still viable. 

Here, we report the Chandra discovery of extended 
X-ray emission associated with the compact steep spec- 
trum (CSS) quasar, 3C 186 (Q0740-h380, z=1.063). The 



Chandra spectrum of the diffuse emission contains 741±40 
counts and a strong (EW~412eV) iron emission line at the 
quasar redshift characteristic of thermal emission. This is 
the first observation of thermal emission associated with a 
CSS quasar at high redshift and gives us a rare opportu- 
nity to study interactions between an expanding CSS radio 
source and the cluster medium. As we show in Section 6.2 
that the pressure of the hot cluster gas is too low, by 2-3 
orders of magnitudes, to confine the radio source. This is 
direct observational evidence that the radio source is not 
thermally confined, and so instead is presumably young so 
that we are observing it at an early stage of its evolution. 

Over the last decade attempts have been made to find 
X-ray clusters associated with radio-loud sources at high 
redshift (e.g. O'Dea 2000, Siemiginowska et al 2003 for 
studies related to GPS and CSS). If the clusters are found 
around large number of radio-loud sources then this could 
be used to place constraints on structure-formations mod- 
els at large redshift. Bremer, Fabian & Crawford (1997) 
describe the model for an onset of a powerful radio source 
in the center of a cooling flow cluster. Such a cluster could 
also confine a compact radio source. The limited capabil- 
ities of the available X-ray telescopes allowed only for a 
few detections of extended X-ray emission around radio 
sources at redshifts z > 0.3 (Hardcastle & Worrall 1999, 
Crawford & Fabian 2003, Worrall et al. 2001). High dy- 
namic range observations are required to detect faint dif- 
fuse emission in the vicinity of a bright powerful source. 
The Chandra X-ray Observatory can resolve spatially dis- 
tinct X-ray emission components in the vicinity of a strong 
X-ray source with ~1 arcsec resolution and a high dy- 
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namic range, as evidenced, for example, by the discovery 
of many resolved quasar X-ray jets (e.g. Schwartz ct al. 
2000, Siemiginowska et al. 2002, Sambruna et al. 2004, 
Marshall et al. 2005). 

The large-scale X-ray emission observed in radio sources 
can result from several processes, confusing detections of 
X-ray clusters at z > 1 (e.g. Celotti & Fabian 2004). 
Radio synchrotron emission implies the presence of a pop- 
ulation of rclativistic particles that produce high-energy 
emission via inverse Compton scattering of the cosmic mi- 
crowave background (CMB) photons. The energy den- 
sity of the CMB increases with redshift as (1 -I- z)'^, so 
the surface brightness of the inverse Compton emission is 
approximately constant with redshift. This is in contrast 
with thermal cluster emission where the surface brightness 
drops with redshift. X-ray spectral and spatial informa- 
tion are key to identifying different emission components 
and measuring diffuse thermal X-ray emission. 

At the highest redshifts Carilli et al. (2002) describe ex- 
tended emission possibly associated with thermal emission 
from shock- heated gas within ~ 150 kpc of the z = 2.156 
radio galaxy PKS 1138-262. Fabian et al. (2003) report 
~100 kpc-scale emission around the nucleus of the redshift 
z =1.786 radio galaxy 3C 294 and give several possible ex- 
planations for the origin of this emission. In both cases the 
poor quality of the X-ray spectrum does not allow confir- 
mation of the thermal nature of the diffuse X-ray emission. 
At redshifts 0.5 < 2 < 1.0 thermal confirmation has been 
possible in some sources, e.g. 3C 220.1 (Worrall et al. 
2001) 

Several deep images of nearby X-ray clusters (Fabian et 
al. 2003, Forman et al. 2003, Nulsen et al. 2004) obtained 
recently with Chandra provide evidence that intermittent 
AGN outbursts with an average power of ~ 10^^ erg sec"-'^ 
supply energy into the cluster preventing its cooling (Mc- 
Namara et al. 2005). These deeply imaged clusters contain 
only relatively low power AGN (~ lO"'" — 10^"^ erg scc"-'^) 
that are interpreted as having been active in the past in 
order to drive the observed cluster morphology. In con- 
trast to these low redshift clusters, the 3C 186 cluster is 
observed during the quasar's active phase while heating 
the cluster medium. 

3C 186 is a very luminous quasar (L;,,,; ~ 10^^ 
erg sec~^). It has a strong big blue bump in the optical-UV 
band and broad optical emission lines (Netzer et al. 1997, 
Simpson & Rawlings 2000, Kuraszkiewicz et al. 2002, 
Evans & Koratkar 2004). It is therefore a typical quasar 
except for its radio properties. The radio morphology 
shows two components separated by 2" and a jet connect- 
ing the core and NW component (Cawthorne et al., 1986). 
Murgia et al. (1999) estimated the age of the CSS source 
to be of the order of ~ 10^ years based on the spectral age 
of the radio source. Our observation provides X-ray mor- 
phology and spectral information for the quasar and an 
associated X-ray cluster. The quasar core is so bright in 
X-rays that any X-ray emission associated with the radio 
components is not spatially resolved. The diffuse X-ray 
cluster emission is detected beyond the quasar core and 

^ http:/ /asc. harvard. cdu/proposcr/POG/index. html 
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the CSS source. 

We describe the Chandra observation, our analysis tech- 
niques and the results for the extended component in 
Sec. 2. Section 3 and 4 present the archival radio and opti- 
cal data. The quasar X-ray spectrum is presented in Sec. 5. 
Section 6 contains the discussion. 

Throughout this paper we use the cosmological pa- 
rameters based on the WMAP measurements (Spergel et 
al. 2003): Hq =71 km scc'^ Mpc'^, fin = 0.27, and 
i^vac ~ 0.73. At z = 1.063, 1" corresponds to ~8.2 kpc. 

2. CHANDRA OBSERVATIONS 

3C 186 was observed for ~ 38 ksec with the Chandra 
Advanced CCD Imaging Spectrometer (ACIS-S, Weisskopf 
et al. 2002) on 2002 May 16 (ObsID 3098). The source 
was located ^ 35" from the default aim-point position 
(to avoid node boundaries) on the ACIS-S backside illu- 
minated chip S3 (Proposer's Observatory Guide (POG)^). 
The 1/8 subarray CCD readout mode of one CCD only was 
used resulting in 0.441 sec frame readout time. The obser- 
vation was made in VFAINT mode with the standard 5x5 
pixel island used to assign the event grades by the pipeline. 
This mode allows for a more efficient way of determining 
the background events and cleaning the background, espe- 
cially at the higher energies. After standard filtering the 
effective exposure time for this observation was 34,398 sec. 
Given the ACIS-S count rate of 0.025 counts s~^ frame"^ 
the pileup fraction was low < 2% (see PIMMS^). 

Figure 1 shows the Chandra ACIS-S image overlayed 
with the quasar core and background regions. Figure 2 
shows the ACIS-S image of 3C 186 adaptively smoothed 
with the CIAO tool CSMOOTH. The X-ray emission is 
more extended than the X-ray emission of a typical radio- 
quiet quasar observed with Chandra. The diffuse emission 
extends up to ^ 120 kpc in radius. Figure ?? show the 
smoothed images in soft (0.5-2 keV) and hard (2-7 keV) 
energy bands. Only the NE quadrant is visible in the hard 
band, while all of the extended emission is visible in the 
soft. Below we describe in detail the imaging and spectral 
analysis of this emission. 

2.1. Imaging Analysis 

The X-ray data analysis was performed in CIAO 
3.2^ with the calibration files from the CALDB 3.0 
data base. Note that the ACIS-S contamination file 
acisD1999-08-13contamN0003.f its was included in our 
analysis; this accounts for the temporal, but not spa- 
tial, variation of the contamination layer on the optical- 
blocking filter of ACIS^. Since the quasar is observed close 
to the aim-point, which is where the spatially-invariant 
contamination model was calibrated, the results will not 
change if the data were re-analyzed using the spatially- 
dependent contamination model. We used Sherpa (Free- 
man et al. 2001) for all spectral and image modeling and 
fitting. 

We ran acis_process_everLts to remove pixel random- 
ization and to obtain the highest resolution image data. 
The X-ray position of the quasar (J2000: 07 44 17.47 H-37 
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53 17.11) agrees with the radio position (Li & Jin 1996) 
to better than 0.1", (which is smaher than Chandra's 90% 
pointing accuracy of 0.6 arcsec, Weisskopf et al. 2003), so 
we have high confidence in the source identification. 

To determine the size of the extended emission we ran 
a ray trace using CHaRT^" and then MARX^^ to create 
a high S/N simulation of a point source. We modeled the 
quasar core as a point source with the energy spectrum 
given by the fitting described in Section 5. We then ex- 
tracted a radial profile from both the Chandra data and 
the simulated point source image assuming annuli sepa- 
rated by 1 arcsec and centered on the quasar. For the 
simulation we added 7% errors to account for uncertainty 
in the raytracc model (Schwartz et al 2000a, Jerius et al 
2004). The PSF was normalized to match the peak surface 
brightness of the core. The resultant profiles are shown 
in Fig 4 and clearly illustrate that the observed emission 
(empty squares) is highly inconsistent with a point source 
(solid line and empty triangles). 

We apply a /3-model to the one-dimensional surface 
brightness profile beyond ^ 3" radius to estimate a core 
radius and /3 parameter for the diffuse emission. The best 
fit model is represented in Figure 5 with (3=0Mtom 

and a 

core radius rc = 5.8"^^ y which corresponds to 47^i[4 kpc. 

Since the core radius found in the one-dimensional anal- 
ysis is small, a two-dimensional fit was made to see what 
influence the quasar emission and any non-sphericity of 
the cluster emission has on the /3-model parameters. A 
model consisting of the ChaRT-generated PSF and a two- 
dimensional (3 model was fitted to the data, using the de- 
fault pixel size of 0.492". The Cash statistic (Cash 1979) 
was used, since the number of counts per pixel was low 
outside the core, and the cluster center, cUipticity, and 
position angle were allowed to vary as well as the nor- 
malization, core radius, and /3 parameter. The best-fit 
cluster location differs from that of the quasar by 0.2", 
which is within the one-sigma error circle of the position 
(0.3"). The best- fit model is elliptical, with an cUipticity 
of 0.24tp;Q7 and position angle of 47 ± 10 degrees, but the 
core radius and /3 values are similar to the one-dimensional 
results, with r^ = 5.5^^]" (45j;Jg kpc) and f3 = 0.581°;°^. 

2.2. Spectral Analysis of the Extended emission 

We extracted the energy spectrum of the extended emis- 
sion from an annulus of radii 2.7" and 15" centered on the 
quasar. There are 1189±34 total counts and 741.4±40.4 
net source counts in this region in the full Chandra en- 
ergy band. The background spectrum was taken from the 
annulus of radii 20" and 30" . Because the background in- 
creases at low and high energies we modeled only the spec- 
trum within a 0.3-7 keV energy range. The total number 
of counts in this energy range was 876±31 with 691.0±32.5 
net counts. In all spectral modeling we simultaneously fit 
background and source data applying statistics with 
the weighting described by Primini et al. (1994) as imple- 
mented in Sherpa. All errors quoted below are 90% errors 
for a single parameter calculated with the pro j ection rou- 
tine in Sherpa. Table 1 lists the applied models and the 
best fit values. 

http://cxc. harvard . edu/chart / 

1 1 http: / /space, mit . cdu/CXC /MARX / 

12 Stark et al 1992 



We first fitted the spectrum of the diffuse emission with 
an absorbed power-law model assuming the equivalent col- 
umn of hydrogen in the Galaxy of 5.68x10^" atoms cm~^ 
(COLDEN^^). We then tested for excess absorption. The 
fitted absorbing column is in agreement with the Galac- 
tic value (see Table 1). Modehng the data with a power 
law indicated an excess at the Fe-line energy character- 
istic of thermal emission at the quasar redshift. We 
added a gaussian line to the model and obtained the line 
equivalent width of EW=412 eV at the observed energy 
Eo6s=3.18±0.07 keV (90% errors) corresponding to the 
rest- frame energy of 6.56±0.14 keV. Note that this line 
energy indicated that the line is emitted by the hot ion- 
ized plasma. 

We next applied the RAYMOND and MEKAL XSPEC plasma 
models in Sherpa and obtained gas temperatures 4.4- 
5.2 keV for two choices of abundance, solar and 0.3 solar, 
using the abundances of Anders & Grevesse (1989). Fig- 
ure 6 shows the best-fit thermal model together with the 
residuals. The observed unabsorbed flux assuming these 
models is of order 6.2±0.3 x lO^^^'erg sec~^ cm~^ for the 
0.5-2 keV energy range and 5±0.7 x 10~"'^^erg sec~^ cm~^ 
for the 2-10 keV energy range (the flux errors are based 
on 5% and 14% counts uncertainties in each band respec- 
tively) . 

The 0.5-2 kcV rest frame X-ray luminosity of the cluster 
emission is equal to ~ 3 x lO''^ erg sec~^ (K-corrected). 
This is only the luminosity calculated based on the emis- 
sion in the adopted annulus extending to about 122 kpc. 
This is about 52% of the total cluster luminosity as esti- 
mated from the radial profile fitting described in the pre- 
vious section. Thus the total cluster luminosity is approx- 
imately ~ 6 x 10^^ erg scc"-'^. 

To investigate a possible non-thermal contribution to 
the spectrum due to inverse Compton scattering of CMB 
photons we fixed the MEKAL model at the best value and 
added a power law component with F = 1.7 to the model. 
The 3(7 upper limit for the contribution of this component 
at 1 keV is then 3.4xl0~^ photons cm~^ sec~^, i.e. less 
than 12% of the extended emission can be non-thermal. 

To investigate possible temperature variations in the ra- 
dial direction we extracted the spectra from two annuli: 
the inner one spans 2.7"-7.8" and the outer one spans 7.8"- 
15". Table 2 summarizes the results. The hardness radio 
indicates that there are more soft counts in the inner re- 
gion than in the outer one. The RAYMOND model fit to 
these spectra indicate a slight temperature decrease, by 
AkT ~ 0.3 keV, towards the inner radii. However the 
errorbars are larger than any apparent deviations. The 
available data do not allow for fitting of cooling flow mod- 
els (see Sec. 6.1 for cooling-fiow discussion). 

To further investigate possible temperature variations 
of the extended emission in the azimuthal direction we 
divided the 2.7" to 15" annulus into four sectors and ex- 
tracted counts and spectra from these sectors. The sec- 
tors, illustrated in Fig. 7a, were chosen to follow the non- 
symmetrical shape of the extended emission as apparent in 
the smoothed image (Fig. 2 and Fig. 3). We calculated the 
surface brightness profile in each sector which confirmed 
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that the emission in the NE-SW direction (sectors 1 and 
3) is stronger (by a factor of 1.25-1.5) than in NW-SE di- 
rection (sectors 4 and 2). The surface brightness vahies 
are indicated in Figure 7b. 

The X-ray emission properties of the sectors are pre- 
sented in Table 2. The hardness ratios indicate possible 
spectral differences between sectors. We fit the Chandra 
spectrum of each sector with a thermal emission model. 
There is a slight temperature variation (with a minimum 
of 3.7 kcV and a maximum of 4.9 kcV) between sectors 
in the best-fit values, however the error bars arc approxi- 
mately 1.5-2 keV, so that the fitted variations are all within 
90% errors. 

In summary, the diffuse emission is non-symmetric and 
elongated in the NE-SW direction with the harder emis- 
sion towards NE (see Fig. 3). Note that this structure is 
orthogonal to the radio emission described in the next sec- 
tion. 

3. RADIO OBSERVATIONS 

We searched for radio emission corresponding to the dif- 
fuse X-ray emission revealed in the Chandra image by re- 
analyzing VLA 1.5 GHz data published by van Breugel 
et al.. (1992). The multi-configuration (A and B) data 
were obtained in 1987 and amounted to a total integration 
time of about 80 minutes split almost equally between the 
two configurations. At the resolution of our image (Fig- 
ure 8a), the radio source is dominated by a 1.8" double 
source aligned at a position angle of about -37 deg. Our 
self-calibrated dataset does not show extended radio emis- 
sion on the scale of the X-ray emission detected above a 
3(T rms noise of 0.75 mJy/beam in the naturally weighted 
image; the off-source (>15") rms achieved in the image is 
about a factor of two smaller. 

We also analyzed a 410 second VLA 15 GHz A configu- 
ration dataset of 3C 186 obtained on 13 Dec 1992 (program 
AL280) in order to examine the radio jet more clearly. The 
image shows similar features to previous high resolution 
maps (Cawthorne et al.. 1986; Spencer et al. 1991), i.e., a 
one-sided jet to the northwest connecting a flat-spectrum 
core to a diffuse radio lobe, and a bright radio lobe in the 
southeast direction (sec Fig. 8b). 

4. HUBBLE SPACE TELESCOPE DATA AND OPTICAL 
EMISSION 

A cosmic-ray rejected WFPC2 Associations^'^ image 
with a total integration time of 8,000 seconds was down- 
loaded from the HST archive. The F675W image, obtained 
in Program 6491, was centered on one of the lower resolu- 
tion wide field (WF) chips. We modeled the central source 
with elliptical isophotes utilizing the ELLIPSE task in the 
STSDAS^^ package and subtracted it from the image in 
order to show the nearest lying objects more clearly. The 
resultant image (Fig. 9) shows several sources within the 
X-ray cluster emission indicated by a 15" circle. Unfortu- 
nately only a single band image of 3C 186 was taken with 
the WFPC2 camera and we cannot identify the colors of 
these sources. 

5. QUASAR X-RAY EMISSION 

^■^ http:/ /archive. stsci.edu/hst/manual 

http:/ /www. stsci.edu/resources/software_hardware/stsdas 



The quasar core emission dominates the overall X-ray 
emission in the vicinity of 3C 186. We define the quasar 
emission region as a circle with 1.75" radius and assume 
background emission from an annulus with radii 20"and 
30" (see Fig.l). Based on the PSF modeling we esti- 
mate that ~98% of the point-source counts are included 
in this source region. The extracted quasar spectrum con- 
tains 1968.7±44.3 net counts (1905.4 net counts in energy 
range between 0.3-7 keV). We model the spectrum with 
an absorbed power law. The best-fit power-law model has 
a photon index P = 2.01±0.07 and a 2-10 keV flux of 
1.7x10"^"^ ergs cm^^ scc^^ which corresponds to a quasar 
X-ray luminosity Lx(2-10 keV) ~ 1.2xl0^^erg sec~^ and 
Lx(0.5-2 keV)=l.lxlO'*5erg scc'^ (unabsorbed and K- 
corrected luminosity). We do not detect any significant 
neutral absorbing column intrinsic to the quasar with the 
3(7 upper limit to the equivalent column of Hydrogen of 
< 9.0 X lO^^atoms cm'^. 

The power-law fit to the data leaves some residuals at 
~ 3 keV indicating a possible emission line at this energy. 
Wc added the Gaussian line component to the model and 
obtained the best fit location for the narrow emission line 
(FWHM< 0.23 keV) at Eobs=i-07toTi keV correspond- 
ing to Erest=6.33±0.06 keV. The line equivalent width is 
equal to EW=162 cV and it can be identified with Fe- 
Ka emission. The best-fit model and residuals are shown 
in Fig. 10. (Note that the residuals between 1 and 2 keV 
are due to calibration uncertainties.) We can estimate 
the flux contribution to this line from the extended ther- 
mal emission by extrapolating the radial profile of the ex- 
tended emission into the central circular region assumed 
for the quasar emission. The contribution from the ther- 
mal cluster emission to the quasar spectrum is of order 
10%. Given the line fiux in both components we esti- 
mate that only about ^5% of the line emission can come 
from the thermal cluster emission, thus the line is dom- 
inated by the nuclear emission. Note that the energy, 
Erest=6.33±0.06 keV, of the detected Fe-line in the quasar 
spectrum is in agreement with being emitted by the neu- 
tral medium, while the energy of the Fe-line emitted by 
the cluster gas, Erest=6.56±0.14 keV, indicates an emis- 
sion from ionized gas. 



6. DISCUSSION 

6.1. X-ray Cluster 

Our Chandra observation reveals X-ray cluster emission 
at the redshift of the quasar 3C 186. The X-ray properties 
of the cluster are summarized in Table 3. We compare 
the cluster temperature and its luminosity with results for 
the other clusters at high redshift using the MEKAL model 
with the abundance set to 0.3 as in Vikhlinin et al.(2002). 
The cluster temperature of ^ 5.2t_Q'Q keV and the total 
X-ray luminosity of Lx(0.5 — 2 keV) 6 x 10^^ erg sec~^ 
agree with the temperature-luminosity relation typically 
observed in high redshift (z > 0.7) clusters (e.g. Vikhlinin 
et al. 2002, Lumb et al. 2004). 

We can estimate physical properties of the cluster us- 
ing standard formulae (Donahue et al. 2003; Worrall & 
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Birkinshaw, 2004). Approximating the X-ray diffuse emis- 
sion as spherical, we calculate the cluster central electron 
density (tiq) to be approximately 0.044±0.006 cm~'^ (er- 
rors are only due to the uncertainty in normalization) for 
the best-fit ID beta model parameters: (3=0.64, a core ra- 
dius of 9c = 5.8" (Sec. 2.1), a gas temperature of 5 keV and 
the spectral normalization of 3.5(±0.4) x 10^* based on the 
MEKAL model fit to the spectrum of the diffuse emission 
(assuming the annulus between 2.7"and 15"radii). These 
parameters imply that the mass of the gas enclosed within 
1 Mpc radius of the isothermal sphere is ~2.2x IO^-^Mq. 
The total gravitational mass enclosed within 1 Mpc is 
~2.6xlO^^M0 assuming an isothermal-sphere model. The 
cluster gravitational mass is comparable to the mass of the 
other high redshift clusters recently measured with Chan- 
dra and XMM-Newton (Donahue et al. 2003, WorraU & 
Birkinshaw 2003, Vikhhnin et al 2002). The gas fraction, 
i.e. ~ 10%, broadly agrees with the gas fraction usually 
found in high redshift (z>0.7) clusters. 

We found a relatively small core radius, ~ 50 kpc, for 
the 3C 186 cluster while most redshift z > 0.7 clusters 
have typical core radii larger than 100 kpc. If gas of the 
same X-ray luminosity as found within a radius of 123 kpc 
of 3C 186 were distributed with a beta model of more typ- 
ical core radius (150 kpc), the gas mass contained within 
a radius of 1 Mpc would increase by a factor of about 2. 

Is there a cooling flow in this cluster? We estimate the 
gas mass enclosed within the core radius of ^ 50 kpc to be 
of the order of 2.2xlO"'^^M0. Given the gas central density 
and the temperature of 5 keV the cooling time for this core 
is ~ 1.6 X 10^ years and without a heat source there would 
be a cooling flow with the cooling-flow rate ^ 50 M0 yr~^ 
(Fabian & Nulscn 1977, Fabian 1994). 

Our image analysis indicates that the cluster X-ray 
emission follows an elliptical distribution (see Section 2.1). 
The angle between the major axes of the two-dimensional 
ellipse and the radio jet axes is ~ 84tg degrees. This 
type of X-ray vs. radio morphology is often seen in lower 
redshift clusters and it is interpreted as due to interac- 
tions between the radio plasma and the cluster medium 
(as for example in Hydra cluster in Nulsen et al 2002). 
The current radio data do not show any radio emission 
on scales similar to the observed X-ray cluster emission. 
However, the X-ray morphology suggests that there could 
be old plasma there possibly associated with previous ac- 
tivity of the quasar. The detection of a relic at low radio 
frequency observation would allow studies of the evolution 
timcscalcs and the feedback between the quasar activity 
and the cluster medium. 

6.2. Optical Environment of 3C 186 

Optical observations of the quasar 3C 186 indicated that 
the source is located in a rich galaxy environment. Sanchez 
and Gonzalez-Serrano (2002) studied an over-density and 
clustering of galaxies in the optical field of several high 
redshift radio sources. They include a K-band image of 
the 3C 186 field and indicate that the surface density of 
the galaxies in a possible cluster peaks to the NE about 
50" from the quasar. While the reported location of the 
peak density is located just outside the Chandra FOV we 
do not detect any significant X-ray emission towards the 
optical peak away from the quasar. In fact the quasar 



seems to be centrally located with respect to the X-ray 
diffuse emission. 

Recent optical study of the cluster environments of 
radio-loud quasars at 0.6 <z< 1.1 by Barr et al (2003) 
shows that these quasars do not reside in the center of 
the galaxy distributions. Because the X-ray emitting gas 
traces the cluster's gravitational potential, the X-ray ob- 
servations can confirm whether the quasar is located in the 
center of the potential well of the cluster. The Chandra 
observations 3C 186 indicate that the quasar is located at 
the center of the X-ray emission and offset from the center 
of the galaxy distribution. 

6.3. Lack of Confinement of the CSS source 

The diffuse emission on a ~120 kpc scale discovered in 
this Chandra observation indicates the presence of a hot 
Intercluster medium (ICM) surrounding this powerful CSS 
quasar. The X-ray spectrum of the cluster is dominated by 
a thermal component with a strong Fe-line at the quasar 
redshift. In one scenario the small size of CSS radio sources 
is associated with a dense environment that prevents the 
expansion of the source and confines the radio lobes to the 
size of a host galaxy (Wilkinson et al. 1981, van Breugel 
et al. 1984, O'Dea et al. 1991). Is this cluster then re- 
sponsible for confining the CSS source? 

Based on the cluster central density and temperature, 
we estimate a central thermal pressure of'^5xl0~^^ dyn 
cm~^. If this pressure is higher than the pressure within 
the expanding radio components of the CSS source then 
the cluster gas may be responsible for confining the radio 
source and its small size. 

The radio source is dominated by a double at 1.5 GHz, 
which straddles a central flat-spectrum radio core as seen 
in the higher resolution 15 GHz map (Figure 8). Us- 
ing DIFMAP's MODELFIT (Shephard, Pearson & Tay- 
lor 1994) routine, the 1.5 GHz data are best-fit with two 

0. 59 Jy elliptical Gaussians with dimensions of 0.7" x 0.3" , 
and a separation of 1.8". Taking the spectral index to be 

1, and approximating each component as a homogeneous 
spheroid, we estimate the minimum pressure in each radio 
component to be ~10^* dyn cm~^ (see also Murgia et al 
1999). Thus the radio source is highly overpressured by 
about 2-3 orders of magnitude with respect to the thermal 
cluster medium. 

The new X-ray detection gives direct observational ev- 
idence that the radio source is not thermally confined as 
posited in the "frustrated" scenario for CSS sources. In- 
stead, at least in this CSS quasar, it appears that the radio 
source may indeed be young (Rcadhead & Hewish 1976; 
Phillips & Mutel 1982, Carvalho 1985) and that we are 
observing an early stage of radio source evolution. 

The jet can also interact and be stopped by a clumpy 
cold medium of the host galaxy, as described by Carvalho 
et al 1998, De Young 1991 or Jeyakumar et al 2005. The 
absorption due to this cold medium should be detectable in 
the quasar X-ray spectrum. We can use the 3(t upper limit 
of 9x10^" cm^^ on the total absorbing column density in- 
trinsic to the quasar (see Sec. 5 and Table 4) to estimate 
the total size of the clumpy medium. The frustrated jet 
models require clouds with densities of 1-30 cm~'^. The 
absorption limit gives the size of the clumpy medium of 
order 10-100 pc compared with the 16 kpc diameter of the 
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radio source. Any such region cannot significantly limit 
the expansion of the radio source and frustrate the jet 
(see also discussion in Guainazzi et al 2004). 

Our X-ray observation of 3C 186 indicates that the CSS 
radio source is not confined, but it is at its early stage of 
the evolution into a large scale radio source. 

6.4. The CSS and the Cluster Heating 

The CSS radio components are overpressured with re- 
spect to the thermal cluster gas. Thus the expansion of 
these components into the cluster medium could poten- 
tially heat the center of this cluster. The energy dissi- 
pated into the cluster by the expanding radio components 
has been widely discussed in the context of the low red- 
shift clusters, where there is evidence for the repetitive 
outbursts of an AGN. However, the details of the dissipa- 
tion process are undecided with ion viscosity and "sound" 
waves being possible candidates (Fabian et al. 2005). On 
the other hand the mechanical energy released during the 
shock wave propagation throughout the cluster could be 
transfered into the cluster thermal energy at the location 
of the shock (P.Nulsen private communication). 

Wc can estimate the energy content of the hot clus- 
ter gas assuming a total emitting volume of 2.3xl0^^cm'^ 
(contained by an annulus with 3 and 15" radii, assuming 
spherical geometry) and kT ^ 5 kcV, to be of the order of 
^kTnV ^ 2x10^^ ergs (where n is the average gas parti- 
cle density in the cluster). We can estimate the jet power 
from the relation between the radio luminosity and the jet 
power given by Willott et al (1999, Eq.(12)). Using the 
151 MHz flux density of 5.9x10^^"' erg sec^^ cm^^ Hz^^ 
(Hales et al. 1993) which accounts for the total radio emis- 
sion from the jet and hot spots (the core is already almost 
completely self-absorbed at 1.7 GHz, Spencer et al. 1991) 
the jet kinetic power is of order Ljet ~ 10"^^ erg sec~^. 
If the expanding radio source dissipated the jet's energy, 
Ljet ^ lO^^crg sec~^, into the cluster's central 120 kpc 
region then the heating time would be ~10^ years. We 
can also estimate the amount of mechanical work done 
by the jet and radio components during the expansion 
to the current radio size (2" x 0.3" ~ 2.3 x lO^^cm^) as 
pdV ~ 10^^ ergs. If the expansion velocity is of the order 
of 0.1c then the radio source has been expanding for about 
5 x 10^ years with an average power of 6x10^^ erg sec~^. 
The estimated jet power is ~ 3 orders of magnitude higher. 

6.5. X-ray emission of the core 

The evolution of radio-source expansion within host 
galaxies and clusters has been considered by Heinz, 
Reynolds and Begelman (1998). They simulated interac- 
tions between a growing radio source and the interstel- 
lar and intergalactic medium. For the highly supersonic 
expansion of the young source a shock forms around the 
expanding source and it heats up the medium to X-ray 
temperatures. As a result a "cocoon" of hot medium sur- 
rounds the radio source. Depending on the density of the 
medium and the strength of the shock a source of the size 
of 16 kpc can emit 10*^ erg sec~^ in the Chandra band. 

The double CSS radio source and the radio jet remain 
unresolved in the Chandra observation of 3C 186. Thus 
the "cocoon" region is located within the unresolved X- 
ray core and thus could be giving a significant contribu- 



tion to the total observed X-ray luminosity of 10 erg s~^. 
The amount of this contribution depends on the physical 
parameters of the expanding radio source and the ISM. A 
weak Fe-line is the only emission feature present in the oth- 
erwise featureless X-ray spectrum, however the observed 
spectral photon index is quite steep for a typical radio 
loud quasar (Elvis et al 1985, Bechtold et al 1994). We 
therefore estimated a possible thermal contribution to the 
quasar luminosity using the RAYMOND plasma model fits to 
the quasar spectrum. We concluded that about ~ 15% of 
the 0.5-2 keV luminosity, e.g. 1.5x10^** erg sec~^ could 
be due to the thermal emission. Of course there are other 
possible contribution to the observed X-ray spectrum from 
the radio jets knots and hot spots. Consistent theoret- 
ical modeling of the expansion and the emission of the 
jet components in the future may help in understanding 
their relative contributions to the X-ray spectrum. This is 
necessary in order to disentangle "true" quasar emission 
related directly to the accretion flow. 

The quasar optical-UV (big blue bump) luminosity of 
5.7xl0^^crg sec~^ (based on measurements in Simpson & 
Rawlings, 2000) is dominated by the typical quasar emis- 
sion related to the accretion onto a supermassive black 
hole. We can therefore estimate the central black hole 
mass and required accretion rate based on that luminos- 
ity. Assuming that the quasar is emitting at the Edding- 
ton luminosity the black hole mass should be of the order 
^ 4.5 X IO^Mq. Based on CIV FWHM measurements 
of Kuraszkiewicz et al (2002) and the Vestergaard (2002) 
scaling relationship for the black hole mass, the estimated 
mass of the black hole is approximately a factor of 10 
higher, ^ 3.2 x IO^Mq. In any case the accretion rate 
required by the observed UV luminosity, and assuming 
10% efficiency of converting the gravitational energy into 
radiation, is equal to ~ 10 Moycar"^. Given the age of 
the radio source of 5 x 10^ years, a total of ~ 5 x lO^M© 
should have been accreted onto the black hole to support 
the current "outburst" . 

Recent studies of formation of galaxies and growth of a 
supermassive black hole suggest that the quasar activity 
is a result of a merger event (e.g. DiMatteo, Springel & 
Hernquist 2005) . In this model a short phase of quasar ac- 
tivity is a direct consequence of the increased fuel supply 
onto a central black hole. Based on the radio aging and 
the current size of 3C 186 its activity started ~ 10^ years 
ago. The X-ray cluster emission of 3C 186, although el- 
liptical in shape docs not show any signatures of a merger 
event responsible for that activity. The cooling time for 
this cluster is 10^ years, so if this cluster is forming then 
the gas flowing into the center could onset the radio ac- 
tivity (Bremer et al. 1997). Such a powerful outburst can 
easily prevent the cooling of the cluster core. Although, 
the exact process of transferring the outburst energy into 
the cluster gas is unknown, the available outburst energy 
exceeds by > 2 orders of magnitude the luminosity of the 
cluster core. The feedback between the jet heating and 
the cooling of the cluster can regulate the growth of the 
central black hole. On the other hand the short timescale 
of the current outburst may also suggest that the intermit- 
tent AGN behaviour could be related to the physics of the 
accretion flow (Sicmiginowska, Czerny & Kostyunin 1996, 
Janiuk et al 2004) instead of the properties of a large scale 
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environment triggering the quasar outburst. This is the 
first cluster observed at the early phase of the quasar ra- 
dio activity giving us a potential to study the early stages 
of interactions between the cluster and AGN. 

6.6. Old halo: relic 

Is the current 3C 186 outburst the first phase of the 
quasar activity in this cluster? The large mass of the cen- 
tral black hole suggests that the supcrmassive black hole 
must have been formed much earlier than 10^ years ago. 
As we show in the previous section only a small fraction 
of the mass has been accreted during the current active 
phase. 

One way to accommodate the observed statistics of 
source sizes in the radio source evolution is to invoke 
intermittent activity (Reynolds & Begelman 1997) with 
average timescales of ^10''^ - 10^ yrs. In this scenario, 
a diffuse radio halo filled with old electrons from previ- 
ous periods of activity should be apparent around CSS 
sources like 3C 186 in deep, low-frequency radio maps. 
Most CSS and GPS sources do not show very large-scale 
(lOO's kpc) extended radio emission in present data (O'Dea 
1998 and references therein). X-ray observations of nearby 
clusters suggest a longer period of the intermittency, e.g. 
10® — 10^ years (e.g. Fabian et al 2003, Forman et al 
2004). For these longer timescales 3C 186 would still be 
at an early phase of the cycle. 

If there was a previous outburst we might be able to 
detect "bubbles" filled with the old radio plasma (seen as 
depressions in the X-ray surface brightness) in deeper X- 
ray observations in the future. The "bubbles" observed 
in nearby clusters indicate locations of radio plasma in- 
jected into the cluster in the past (Churazov et al 2001). 
They are buoyant, moving within the cluster gas and last 
for at least one cycle, as indicated by more than one pair 
of "bubbles" observed in some X-ray clusters (Fabian et 
al. 2002, Belsole et al 2001, Owen et al. 2000). Even if 
the radio source outburst lasted only ~10^ years the radio 
plasma would still be present within buoyant "bubbles" 
loosing its energy via synchrotron radiation. Thus if there 
was a previous outburst we might be able to detect the 
radio relic in low-frequency observations. 

At large redshifts, large-scale radio halos may produce a 
non-thermal X-ray component via inverse Compton scat- 
tering off the CMB. In the case of 3C 186, at z=1.063, 
the X-ray spectral fits do not require such a power-law 
component in addition to the thermal emission. From 
our data (Sec. 2. 2), we estimate that a non-thermal com- 
ponent contributes less than 12% of the total emission. 
F(l keV)-few nJy. 

No large-scale extended radio emission is visible in our 
1.5 GHz VLA image - we estimate that there is less than 
0.3 Jy coming from a possible halo by integrating the 3cr 
rms limit over the extent of the detected X-ray emission 
(~15"). The electrons scattering CMB photons to 1 keV 
(7 ^10'^) radiate synchrotron radio emission at much lower 



frequencies (j/ (B / 1 ^G ) MHz). A deep low-frequency 
radio observation aimed at detecting extended radio emis- 
sion on the scale of the X-ray cluster would put a useful 
lower limit on the cluster magnetic field assuming equipar- 
tition (see CariUi & Taylor 2002). 

7. SUMMARY 

Our main results can be summarize as follows: 

1 We have observed the CSS quasar 3C 186 for 38 ksec 
and detected X-ray cluster emission extending out 
to ~ 120 kpc from the quasar. 

2 The cluster temperature and luminosity follows the 
relationship observed in the other high rcdshift clus- 
ters, its gas mass is relatively small. The low mass 
could be the result of a small core radius measured 
in this cluster in comparison to the other clusters 
at z > 0.7. 

3 The estimated pressure of the cluster gas is 2-3 or- 
der of magnitude lower than the pressure of the ra- 
dio components, thus the cluster gas cannot confine 
the expanding radio source. 

4 Non-thermal and thermal emission associated with 
the radio components cannot be separated spa- 
tialy from the unresolved X-ray emission that is 
measured. We have placed an upper limit of 
1.5x10"** erg sec~* on the contribution from the 
gas that is shock heated as a result of the radio- 
source expansion. Detailed modeling is required to 
estimate the relative contribution of several possible 
components identified in radio observations. 

5 Future low frequency radio observation may provide 
information about the large scale distribution of an 
old electron population from a previous outburst of 
the quasar activity in this source. 
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Table 1 

Spectral Models of the X-ray Diffuse Emission 



Model" 


N/f or z 


Abund r 
Solar 


Temp 
keV 


Norm'' 


X2" 
(DOF) 


PL 




o nc + 0.27 


- 


3.08t«fg 


79.92 (73) 


PL+Gaussian Line 
Ef = 3.18±0.07 
EW=412.7cV 




2.12±° 12 




3.06±° 21 


68.69 (71) 


RAYMOND 




1 


^■°-0.7 


30.2t^j 


73.7 (74) 


RAYMOND 




0.3 


. 2+1.3 
^■^-0.9 


34.8t|;i 


74.0 (74) 


RAYMOND 


1 11+0.04 
^■^^-0.03 


0.3 


5 4+1'' 
'-'■^-1.0 


36.3l^t 


71.1(74) 


MEKAL 




1 


4 6+11 


29.8i^;^ 


71.9 (74) 


mekal 




0.3 


r 2+1.3 
^■''-0.9 


34.6_3 


73.7 (74) 



" All models except the first PL model include equivalent Hydrogen column of 5.68x10'^'^ cm ^ in the Milky Way from COLDEN; ^ Normalization in 
10'""' photons cm"^ scc^"'" at 1 kcV for a power law model; for thermal models the normalization is defined as Norm X 10""/(47r(D(l + z)f) J 
following the definitions of RAYMOND and MEKAL models in XSPEC, abundance table set to Anders & Grevcssc (1989); ■= X following Primini ct al 
in Sherpa; ^ Energy of the emission line in keV. 



Table 2 

Spatial Models of the X-ray Diffuse Emission 



Region 


Counts 
Tot 


Net Counts 
Tot 


Net Counts 
0.3-7keV 


Soft 
0.5-2keV 


Hard 
2-lOkeV 


Soft /Hard 


Entire Annulus 
2.7"-15" 


1189±35 


741.1±40.5 


691.0±32.5 




143.91-;^ 




Inner Annulus 
2.7"-7.8" 


490±22 


392.9±24.2 


379.3±21.4 


317.2t_l°-l 


63.2t-« 


r n + 0.9 
°-^-1.2 


Outer Annulus 
7.8"-15" 


699±26 


348.2±32.4 


311.6±24.5 


247.11-- 


o, Q+16.9 
82.8_jg 3 


o 1+0.4 
■5-1-0.6 


Sector 1 


206±14 


142.1±16.4 


128.1±13.4 




37.6i«i 


9 S+0.7 
^-*'-0.9 


Sector 2 


430±21 


272.2±24.2 


246.1±19.3 


200.9tll° 


Kn K+11.'' 
60.6_i2.2 


3 4+°-'^ 

■5-^-0.7 


Sector 3 


164±13 


111.2±14.7 


107.6±12.3 




11 q+5.5 


10 4+1-1 


Sector 4 


39.5±20 


230.8±23.6 


219.6±18.8 


181.4i;?| 


34 7+10.5 


°- ' -2.2 



Table 3 

Summary of Spatial and Spectral Fits to Diffuse X-ray Emission 



Parameter Property 

/3-model (ID) /3=0.64+^-^^, rc=5.8lf ^" 

/3-model (2D) /3=0.58t°:[J^, rc=5.5tli" 

ellipticity = 0.24tj] !]^, PA=47±10 degrees 

Kbs (Fc-line) 3.18±0.07 kcV 

EW (Fe-line) 412 cV 

Fobs(0.5-2 keV) 6.2 ± 0.3 x 10"" erg see"! cm-^ 

Fobs(2-10 keV) 5.0 ± 0.7 x lO^^^ ergs sec-^ cm-^ 

^nonthermal (1 kcV) <5.4 X IQ-l^ erg seC~l Cm"^ 

Ltot(0.5~2 keV) 6x 10^^ erg sec'^ 



Fluxes are unabsorbcd. Luminosities are K-corrcctcd and in the source frame, (see text). 
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Table 4 
Quasar Models 



Modcr 


N(z,,„)'' 


r 


Flux'' 
[1 kcV] 


[koV] 


EW,i„e 

[cV] 


(DOF)'' 


Power Law 


< 9.0" 


2.01±0.07 


7.52±0.32 


+ 0-07 

3.07-016 




138.2 (137) 


PL+Gaussian Line 




2.03±0.07 


7.48±0.28 


130 


132.1 (135) 



^ model assumes 5.68x10^" atoms em~^ equivalent Hydrogen column in the Milky Way from COLDEN; 
^ absorbing column intrinsic to the quasar at z— 1.063 in units of lO^'^ atoms cm~^; 
flux in units of 10 photons cm^" sec^ keV^ X following formulation of Primini et al. in Sherpa: " 3(T upper limit. 



o 




Fig. 1.— Chandra ACIS-S X-ray (0.3-7 kcV) image of 3C 186. The pixel size is the standard ACIS-S pixel of 0.492". The quasar source 
region assumed for extracting the spectrum is a 1.75"radius circle in the center of the field. The outer annulus with 20" and 30" radii 
illustrates the background region. The arrow in the upper right corner indicates the 10"scale. North is up and East is left. Coordinates are 
J2000. 



Fig. 2. — Adaptively smoothed exposure corrected image (photons energies within 0.3-7 keV range) of the Chandra ACIS-S observation of 
3C 186 (Q0740-I-380) The diffuse emission is detected on > 100 kpc scale, 1"=8.2 kpc. North is up and East is left. Contours represent a 
surface brightness of: (0.046,0.066, 0.13, 0.2, 0.33, 0.46, 0.66, 6.635, 33.175) X lO"'' photons cm-^ arcsec^^ -pj^^ direction of the CCD readout 
is indicated by arrow on the right side. A blue arrow in the upper right corner shows the PA=-37 deg of the 2" radio jet (see Fig 8. 



Fig. 3. — Adaptively smoothed exposure corrected image of the Chandra ACIS-S observation 3C 186 (Q0740+380) in two X-ray bands: 
(1) soft 0.5-2keV on the left and (2) hard 2-7 keV on the right. North is up and East is left. Contours represent a surface brightness of: 
(0.036,0.066, 0.11, 0.3, 0.66, 6.635, 33.175)xl0"'' photons cm'^ arcscc-^. i"=8.2 kpc. 



X-ray cluster at z= 1.063 
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Fig. 4. — The background subtracted surface brightness profile up to 25"distance from the quasar. The source data is indicated by the 
square points. The solid line with triangular points represents the simulated PSF profile. 
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Fig. 5. — Background subtracted surface brightness profile for radii 3" to 30" fit with a beta modeL The data are indicated by the square 



points. The solid line shows the best fit model with parameter /3 = 0.64^Q j!,^ and a core radius of rcore 
illustrates the differences between the data and the model in units of a. 



5.8 



11+2.1 



The bottom panel 



X-ray cluster at z= 1.063 
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Fig. 6. — ACIS-S spectrum of the diffuse emission between 2.7 and f 5 arcsec which we fit with the plasma model. Upper panel shows the 
data indicated by squares with la error-bars and the best-fit model drawn with a solid lino. The bottom panel shows the residual difference 
between model and the data in units of sigma. The scatter is due to calibration uncertainties. 
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Region 



Fig. 7. — Left: Defined sector regions overlayed over an Chandra ACIS-S image of 3C 186. One pixel corresponds to 0.492". North is 
up and East is left Right: Background subtracted surface brightness calculated for each sector in the left panel. The sector numbers are 
indicated on the bottom axis. 




Fig. 8. — Left: VLA 1.5 GHz image of 3C 186 from reprocessing the A- and B-configuration archival datasets published in van Breugel 
et al. (1992). The restoring beam is 1.62" X 1.44" at position angle -42.7 degrees shown at bottom left. The image peak is 565 mjy/bm 
and contour levels begin at 0.5 mjy/beam (2cr) and increase by factors of \/2. Some extended radio emission is apparent, though not at the 
angular scale of the observed extended X-rays. North is up East is left. Right: High resolution (0.15") VLA 15 GHz image of 3C 186 showing 
the core-jet morphology of the source. The image peak is 21.6 mJy/beam, and contours begin at 0.65 mJy/beam increasing by factors of \/2. 



X-ray cluster at z= 1.063 
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Fig. 9. — Left: HST image of the field. The 15" radius circle highlights the region of the diffuse X-ray emission. North is up East is 
left. Right: The host galaxy contribution subtracted from the central image. The ring of diffuse emission between 2-3" is an artifact of the 
imperfect galaxy subtraction. Also, the four diffraction spikes (the cross at the center) are apparent. 
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Fig. 10. — Upper panel shows the ACIS-S spectrum of the z=1.063 quasar Q0740+380 (empty squares) over-plotted with the best-fit power 



law plus Gaussian line model at Eoj, 



: 3.07 



+0.06 



keV (a solid line). The bottom panel shows the difference between model and the data in 



units of sigma. The scatter between 1 keV and 2 keV is due to calibration uncertainties. 
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